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ELECTROHYDRODYNAMIC INDUCTION PUMPING 
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by 


Millard Sherwood Firebaugh 
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Engineering and to the Department of Electrical Engineering 
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for the degree of Master of Science and Naval Engineer. 


ABSTRACT 


An electrohydrodynamic traveling wave induction interaction with 
the bulk of a slightly conducting liquid (conductivity of approximately 
10° nos /meter ) is shown to produce a fluid flow. A gradient in con- 
ductivity normal to the direction of flow is required. In the scheme 
described here the conductivity gradient is provided by taking advantage 
of the temperature dependence of conductivity. A traveling potential 
wave is created on an arrangement of electrodes parallel to the flow 
direction and in contact with the liquid. The resulting fields induce 
pcs in the bulk which lag the traveling potential wave. Therefore, 
a time-average electric traction in the bulk is created, motivating the 
liquid. Expressions for the fields, the time-average traction, and the 
fluid velocity are derived and discussed. Experimental results are shown 
which lend credence to the derived equations. 
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CHAPTER I 
INTRODUCTION 


Induced charges in a slightly conducting liquid can interact with 
traveling wave electric fields to produce a fluid flow. One method of 
producing this kind of electroconvection is to introduce a potential 
wave traveling parallel to an air liquid interface. Induced charges 
in the liquid will relax to the interface and form a traveling wave 
of surface charge which will lag behind the potential wave. An electric 
surface shear will result and the liquid will Parte in bors kind ol 
device the interface is a singular point in a gradient of conductivity. 

An extension of the concept of the surface interaction would be 
GO produce a continuous gradient in conductivity in a liquid. Then 
eliminate the interface by bringing the electrodes carrying the travel- 
ing potential wave into contact with the liquid. Charges induced in 
the liquid will no longer relax to the surface but will remain in the 
bulk of the liquid. A kind of internal wave of induced charge will 
then travel along behind the traveling potential wave and the possibility 
for an electric traction in the bulk of the liquid is created. 

This phenomenon is very dependent on the charge relaxation time 


which is the ratio of permittivity to conductivity. If the charge 


* 
Melcher, J.R., "Traveling Wave Induced Electroconvection" to be published 
in Physics of Fluids. | 
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relaxation time for the liquid is very short, then charge will be 


induced in the Liquid but will relax so quickly that it will not 

lag behind the traveling potential wave. On the other hand, if the 
relaxation time is very long, then very little charge will be in- 

duced in the material and the interaction with the traveling field will 
be negligible. 

EHD induction pumping would be feasible for a wide class of 
liquids whose relaxation times would be compatible with reasonable 
frequencies for which traveling electric potential waves might be 
generated. 

In succeeding chapters, a theory and experiment will be described 
which deal with the phenomenon of electrohydrodynamic pumping of liquids 


in the bulk. The results obtained are somewhat surprising. 


=e = 


CHAPTER 11 


THEORY 


AR TAG TOAUC LON 


The basis of the following theory for electrohydrodynamic induction 
pumping in the bulk of a slightly conducting liquid is. that the movement 
of charge in the liquid may be explained by a pure conduction model. This 
model applied to induction pumping in the bulk results in currents which 
are small enough that the problem may be considered to be only an electric 
problem and magnetic effects may be ignored, The theory as developed 
here leads to a conclusion which at first glance is somewhat startling 
and violates an intuitive conception of how an induction device ought to 


operate. 


B. Electrical Theory 
1. Electric Field Equations 
Considering magnetic effects to be small we may immediately write 


k 
Maxwell's equations in the bulk of the liquid. They are: 


vxE=0 (1) 


Vc eB=q (2) 





* 
Fano, Chu, Adler, Electromagnetic Fields, Energy, and Forces, p.179. 





where q is the free charge, e is the permittivity and E is the electric 


field vector. To satisfy conservation of charge 





= Ò 
Vid + = = 0 (3) 


The free current J results from Ohm's Law conduction which written in 


the coordinate system fixed to the laboratory will be 


J =oE + av (4) 


where o is the electrical conductivity and v is the velocity in the 
fluid flow. Since the Electric Reynolds number given by the formula 
"Re = ev/oa is small for the problem under consideration (‘Re ~ 2.6 x vo) 


the qy term in Eque (4) may be ignored leading to 
J = of (5) 


Equation (1) allows that E may be written as the negative gradient of 


the potential 6. 
E=-Ma (6) 
Combining Equs. (2),(3) and (6) leads to: 


9+ (cV8)+S(V +(e V8) )=0 (7) 


ee 





If both o and e are functions of y in the coordinate system shown 
in Mig. 1, this equation is in general difficult to solve, however, 
there are situations in which it may be easily approached. 

One such case occurs if it is assumed that e is a constant 


throughout the bulk of the liquid and that o may be represented as: 
o=o to (c -= 6) (8) 


where € is the normalized depth in the channel. € = y/a. Co re- 
presents some appropriate depth in the channel where a value of o 
is measured which is in some sense an average value of conductivity 
in the channel. It is also necessary that the values of Oy and To 


be such that it is reasonable to say that o ~ o in che Liquids bun 


Vo = o,/4, Where —— ls tie Uni vector am the y direction. 
ly ly 
Combining equations (8) and (7) with e a constant leads then to: 


2 | 5 = 
a, Y a aa a aos (9) 


This equation will now be solved for the two-dimensional case. 


2. Traveling Wave Solutions 
If the upper surface of the channel is excited with a traveling 


potential wave of the form: 


A j (we - kx) 
V=ReVe (10) 





j (wt - Rx) 
Y 2 Re Vo € 





Figure 1 Coordinate system in the pumping channel 


Ge 





where V is the magmi tude of the applied voltage, wis the angular 
Erequency, K iS the Wave number, and circumflex , indicates a complex 
city then the assumed form of the solution of Equ. (9) is: 
A j(wt = kx) 
$ = Re $ (5)e (11) 


Substituting Equ. (11) into (9) gives: 


¿2 


ajina pean 


ti > 


O x A 
> : Sie oO ie) 
C 7 + JWe OC 


7 








o, 


This equation may be solved by assuming for ® solution of the form 
AeP®, Substituting this solution into Eau. (12 ) and solving for 


the roots gives: 


u O O D 
p=3 - —— E ss waxy (13) 
. 9, + jue o + jwe e 


Since in the experiment performed in Chapter IV, „(ax)“ << (o, (6, - 
2 

jwe)) this square root may be expanded around the first term, At 

the Same time defining new constants ù = 9/5, and S = we/o_ the 


roots may be rewritten: 


en 


| 2 ; 


| 





A 
Now the complete solution for ¢ may be written and the boundary 


condition on the electric field in the channel may be applied to 


| evaluate the constants of integration: 

| A p, 6 Pob 

3 (6) =A€ 7 +A, (15) 
At the perfectly conducting bottom of the channel where £ =0, 
the appropriate boundary condition is that the tangential com- 
ponent of electric field must be zero. Therefore, -08/a& = O 
which says that A, = Ay. At the top of the channel where 6 =1 
the condition is that © = V. This condition says that: 

Y, 


AFD D 
el-e? 


The final form of is: 
DS 
Re JANG ee ) 


j(wt - kx) 
6 = 6 


(16) 


Since E, = - 66/a and Le 38/43í. , writing E, and E, in the 
form; 


j(wt - kx) „ d(wre - kx) 
; BE =Re ES 





leads to: 





Se: De 
~ dk V(@ i 
A en 
1 2 
NE E 
A P E Pp E 
A V o (PE l = Poe = ) 
E = 19%) 
y 4 Poo ë 
(E re, 


Equs. (17) and (18) complete the solution of the electric field 
problem for the fields in the bulk of the liquid. 
zz Via Theory 


l. Force Equation 


The force equation will be written for a two-dimensional 
flow in a re-entrant channel under the influence of an internal 
Ale curic shear force which is constant in the x-direction because 
only the time-average value of the electric stress is considered. 
The liquid is incompressible and the flow is steady. 


The general force eugation is: 


een nS ag) EPEN T - pg = (19) 
at : 





P is the pressure, v the flow Melocm sp density, sand Zei, 


the gravitational acceleration. ce is the total stress tensor 





which includes both the electric stress and the viscous stress. 
Under the condition described above, the term p(dv/dt + 


(v . Viv ) = 0 and Equ. (19) may be re-written as: 


er oT 











SEN A (20) 
OK cx Oy 
3 Be OF 
Of = yx + eG m ee 
dy ok oy 
E V r E 
The total stress T.. =T.. +T.. . The electric stress T., = 
lJ LJ LJ LJ 
X 
a < = A 
a 5 si e BB. and the viscous stress has Al (dv, /axj + 


xy 
ev /8x, ), where p may be a function of space. Substituting these 
J 


expressions into Equ. (20) leads to: 








B A E oy (EE ) dv dv 
BA e co ( x x) + u x +e Xx y + al x +) (22 ) 
Ox 2 Ox OK dy dy vx 


In the re-entrant channel there can be no gradient in the x 
direction in pressure. All of the derivatives with respect to x 
in Equ. (2°) are zero because if we consider only time average 
eerie Stresses, chen all of the stresses in the liquid are 


constant in the x direction. This leaves only the last two terms 


ico dson and Melcher, Fields, Forces and Motions, Ch. 5, p. 35. 
** Chandrasekhar, Hydrodynamic and Hydromagnetic Stability, p.le. 


-] OQ 





of Equ. (22) which is the same as saying: 


oT 
a2 = 0 (23) 
oy 





which may be immediately integrated to give: 


T  =C (24) 


where C- is the arbitrary constant of integration. 

Equ. (2L ) will not be treated since it is only concerned with 
the balance of electric forces and gravitational forces in the liquid. 
This equation concerns the internal instabilities in the flow but 
me-oeaQu eiteet the flow in the x direction, at least in the approxi- 
mavlons ot this theory. 

Peene Era ge Or bie electric Stress. which hes neretolore 


only been alluded to, will now be introduced, resulting in the follow- 


ing equation: 


= 


ee ee = 


Q7 
SG 
PS 
9 
er) 
50) 
O 
bi > 
rs 
td > 
* 





The pulsating double frequency term which results when multiplying 


two sinusoids has been ignored. 


-lLle 





2. Computation of Time-Average Electric Stress 
Certain algebraic tricks will be used in this section to perform 
Some computations which are essentially quite tedious. The roots Pp 


and p, given by Equ. (14) may be rewritten as: 


P= | a + | (26) 


+ ja, = 17 /(1 +38) 
EST je) iat wel + jS)/ 


From Equs. (17) and (18) it is possible to write that: 


where a = 


and b 


x x 
rary PS Pob * DS y PS 
Mee Feu Vy (6 = AO (27) 
Re RE = - qqa aan 
0 Bo Pe 
Cc Vee ES) 


where the asterisk indicates the operation of taking the complex 


conjugate. 


p P 
Consider now a term such as € e € = Re-writing after substi- 


tution for the p's from Equ. (26) leads to: 


p P -a/ 
E 


Now sirh b/2 = sirh (b /2 + jp /2) which may be expanded by double angle 
r 1 


formulas leading finally to: 


Da 





Pp -a/2 b b 
iL 2 : : b 
Coa 6 ° = eC (sinh Y cos + j cosh 


= 
2 


>| 


By manipulative tricks of this sort with which the reader may not 
wish to be bored it is possible to reduce Equ. (27), not, however, 
ee Out some effort, to the following form; 


2 -a (6-1) 


kV e . a . = 
se ae lb sin bie b sinh dS + a, (cosh dG cos », 0) (28) 


A A 
Re EE = 
x Vy 
24 (cosh b_- cos b,) 
r í 


The values of the a's and b's may now be computed by combining Equ. (14) 
with Equ. (26). In the experiment described in Chapter IV, it turns 

out that n / (1 + só) >> (ax)*/n. Ignoring the term on the right side 
of this inequality in the expressions for the p's results in the follow- 


ing values for the a's and the b's: 


Em n/a 43°) and T b, ee) ee só) 


It should be noted that taking the above values for the a's and b's is 
completely equivalent to assuming that one solution for the potential 


ö(C) is that @ is a constant. 
3. Velocity Equation 


Equ. (25) may be integrated to determine the velocity. The values 
of the a's and b's are substituted into Equ. (28) and that result substi- 


tuted into Equ. (25), giving: 


232 





A 





— 





n (1-6) NC 
f 














1+S* > E 2 
2 Sl i 
kv" NS e + sin (———-) + cos GA, =e (+S) } Gad 
o `’ S ee Bo 1 
pp a aka a Naane aaa E (29) 
2 N on 
(1487) u(t) (cosh ( > - cos ( 5) ) ulG) 
1+5 1+S 


For purposes of integration and evaluation the boundary conditions on 


the flow at Ç = O and € = 1, Equ. (29) may be written as: 





dv u nd 
ee 1 (30) 
ag me) se uc) 
where: 
ekV TS 
D = _ A MMMM. (31) 


(148°) (cosh E) - cos —) ) 


and: 











A). e f == sin 5) + cos ( ao } (32) 


StS 





Integrating Equ. (30) gives: 





E E 
- = VEN as Po 
(6) » (6) 
O O 


mie boundary conditions for the flow are that at C = 0 ama € = 1, 
the bottom and top of the channel, Voge O a i Chapter TIL it iS 
shown that u(6) = p (7 where n is a constant and this gives (when 


substituted in the above integrals): 


c G 
o. (E) = - (2 P(C) ¢* ag + \ cial 
A Fo 
O O O 


which automatically satisfies the boundary condition at Ç = 0. Ca 


may be evaluated by satisfying the boundary condition at £ = 1. This 


leads to: 

a |, 
MA NENG aci EN (eye al (33) 
LO) ee 


mae velocity equation has several properties of interest, some are 
very interesting, in fact. First, it may be noted that the velocity 


2 
is proportional tovo which appears in the constant D, and that when 


A 





Vy = O the velocity will be O. This is the usual result in electric 
field problems of this sort. Second in the limit as S — O which is 
the limit for O frequency £(6), goes to some constant value and D— 0 
E@mumere can be no motion for O frequency.” The Timi foru = S miy 
not be taken since it violates the conditions for which the equation 
was derived. 

Surely the most interesting result of Equ. (33) concerns the de- 
pendence of vy on TN. The sign of TM is determined by whether or not 
the conductivity gradient in the material from C = O0 to 5 = 1 is positive 
or negative. A positive gradient corresponds to 71 > O and a negative 
gradient corresponds to N < 0. In the expression for f(€), Equ. (32), 
anche Sign of Mus changes ‚the only ua changes its sign is the 
sin term. When the expotential term € (32) dominates the expression 
in brackets in Equ. (32) as it will most of the time,a change in the 
sign of Ç will not change the sign of f(G). Therefore, the expression 
in brackets in Equ. (33) will not change sign . However, the constant 
D will change sign . That is, with a change in’ Sign of the conductivity 
gradient, the pump may be made to move the fluid in a direction opposite 
to that of the traveling potential wave. This amounts to an inverse 
induction device. This direction reversal depends on the sign of the 


puree inaucea in the diguid. Ii tha conductivity gradient is positive, 


mae ts une Fluid is more Conducting at the top than at the bottom, then 


2]6= 





the charge induced in the liquid is positive. (Here we assume that E 
is negative,) The field is slightly ahead of the induced charge because 
of the charge relaxation in the liquid and the resulting electric force 
On the charges is opposite to the direction of travel of the field. I, 
on the other hand, the conductivity gradient is negative, then the charges 
induced are negative and the field tends to pull the charges in the direction 
© the Tield, dragging the fluid along with the charge. 

The dependence of the sign of the charges on the sign of the conductivity 
gradient may be easily shown for the static case. Consider a liquid 
between two parallel plates with a voltage +V | applied across the plates 
from y = 0 to y = d, where " is the direction normal to the plates and a 
gradient in o exists in the liquid such that o = o_ +(0, /a)(y - yo) In 
this case combining Equ. (2), (3) and (5) and noting that dq/dt = O for 


the static case it follows that: 


O 
E 
ge Eo. a oy | (34) 
on d 


where E, is a constant determined LEST. In this static case a positive 
gradient in o corresponds to a positive charge, and changing the sign of 
of JUN changes the sign of the charge. 

Even though at first glance a potential wave traveling in the +x 
direction inducing a fluid flow in the - direction may seem startling 


the theory predicts such a phenomenon and the static case would seem to 


ee 





Support this concept. 
In Appendix A a dynamic problem with a singularity in con- 
ductivity gradient is done exactly and supports the conclusion of 


This theory. 


aos 





CHAPTER IIL 


FLUID PROPERTIES 


A. VrodmeLIon 


The liquid used in the electrohydrodynamic pump was AROCLOR 1232, 
manufactured by the Monsanto Chemical Compeny. The properties of the 
liquid which are of interest in the EHD pump experiment are the electri- 
cal conductivity o, the dielectric constant, and the absolute viscosity 
u. The value of the dielectric constant for AROCLOR was taken from the 
manufacturers published eee The conductivity and the viscosity were 


determined by experiments. 


ibae Conductivity 


Measurements of the conductivity of AROCLOR are not simple to perform 
in terms of repeatibility of data. The apparatus used in the experiments 
was relatively straight-forward, however. The conductivity was measured 
in a so called conductivity cell shown in Fig. 2. The cell consisted of 
a lower plate of aluminum, an upper plate of brass and pyrex walls. To 
measure the conductivity it was only necessary to place a sample of AROCLOR 


in the apparatus so that it completely filled the region between the plates. 


X 
AROCLOR Platicizers, Monsanto Chemical Company, Technical Bulletin, 
No. PL = 6 3 Ds 


2192 





A voltage was applied to the top plate. A 5 megohm resistor was 
connected between the bottom plate and ground and the voltage across 
this resistance was measured with a very high impedance (100 megohm ) 
probe as shown in Fig. 3. Calling the value of the 5 megohm resistor 
Rw the applied voltage V, the voltage measured across the resistor Na 
the separation of the plates d, and the area of the plates A, the con- 
ductivity o was computed from the formula: 


a V 
Y 





In order to determine the temperature dependence of the conductivity 
the apparatus described above was placed in a bath so that the level of 
the bath was above the level of the liquid inside the plates. For 
measurements above room temperature, this bath was mineral oil, heated 
with immersion resistance heaters. For measurements below room temperature, 
the bath was alcohol with dry ice. A thermometer was mounted in the 
apparatus so that the bulb was immersed in the AROCLOR, When the apparatus 
was placed in the cold bath frost tended to form on the sides of the 
pyrex shorting the top plate to the bottom plate. Consequently, the top 
was separated from the pyrex by blocks of plexiglass. This made the separa- 
tion between the plates Large enough for fringing fields to increase the 
effective area of the plates. This effect was corrected for by making 


measurements at room temperature with both the large and small separation, 


-20- 





Thermometer 
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A 
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N E : 
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2h eS eee IT) AS 
Aluminum bottom 
igure 2 Conductivity cell for measuring conductivity of 


AROCLOR as a function of temperature. 


Cell 
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A Conauctivity cell eircuit. 
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comparing the resulting values for o and calculating a correction 
factor. 

The results of the above experiment are shown in Fig. 4. 
aaecse measurements were done at D.C. with a voltage applied of 
6.0 kilovolts. The plate separation d was 2.15 em, A was 33.6 sq. 
em. and the correction factor for fringing was k = 0.815 in the 
Wa” Corrected * i measured’ 

The nature of the actual mechanism for which electrical con- 
ductivity is a manifestation seems to be quite complicated in AROCLOR. 
By placing a point source of light in a position so that the Light 
Shines through the liquid and falls on a white screen, motion of 
the liquid may be easily observed. In an experiment performed with 
the conductivity cell described above with a temperature gradient 
applied across the cell, the motion of the liquid was observed as 
a function of applied voltage. In this experiment the temperature 
of the AROCLOR on the bottom was about y and on the top about 
25 C.: The plate separation was about 2 cm. for a gradient in 
temperature of 25°C /em. With an applied voltage of less than 5 


kilovolts at 1.5 cps. very little motion was observed. From > kilo- 


volts to 7 kilovolts some local swirling of the liquid was noticed 
although it was not violent. At about 7 kilovolts a vaguely defined 


internal surface formed which undulated slowly. Above about 8 


22 





conductivity in mhos/meter 


20 +x 199 


CONDUCTIVITY OF AROCLOR 


VS 
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Lha iji E A A A 
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Figure 4 Results of measurement of conductivity of AROCLOR as a 
function of temperature for DC voltage of 6 kilovolts. 
Lower scale indicates conductivity as a fynetion of nor- 
malized depth in pumping channel. 


ee 





kilovolts this surface was completely eradicated by moderately 
violent swirling. Near 11 kilovolts the pattern of motion changes 
radically to very violent swirling. Under the conditions of this 
experiment gravity would tend to be stabilizing the liquid since the 
colder, denser material is on the bottom and the warmer,less dense 
material is on the top. However, the electric field is trying to 
are che material. Presumably, the electric field@torces over- 
balance the gravitational force at some level of voltage and the 
fluid becomes unstable. This picture is complicated by the varying 
charge density within the material. 

The result of the internal instabilities is to complicate under- 
standing of the mechanism of conduction in the liquid. Mechanical 
motion of the liquid caused by the presence of the field, so-calied 
electroconvection, may be very active in the process which results in 
current flow at the terminals of the conductivity cell. 

In spite of the complexity of the mechanism for conduction in the 
AROCLOR, the data presented in Fig. 4 is assumed to represent the con- 
ductivity in the liquid as a function of temperature. The theory pre- 
sented in Chapter IIT requires that o be represented in the liquid as 
0, T S (€ ~ er As will be explained in Chapter IV, the tempera- 
ture in the liquid in the pumping channel was assumed to be a linear 


function of the normalized distance € with the cold bottom corresponding 


SDN 





to 6 = O and hot top corresponding to 6 = 1. 


Several approaches might be envisioned for approximating the 
data of Fig. 4 as o = O (G~ Cole The first would be to say 


that To is the mean value of o and that T is the mean value of 


do/d6. In that case, co = 10.5 x ora mho/meter and c a 


1 = 129 x 170% 


mho/meter. 6 would then be 0.89. A second possibility would be to 


say that A and co. should be evaluated at G = 0.5, that is, at the 


als 
mid-depth of the channel. In that case o = 2.8 x or mho/meter 
and On Is No mho/meter. Of course, 6 = 0.5. A third 
possibility would be to call on the average value of g and f, the 


TO 


average value of do/d¢. In that event, o, = 4.0 x 10 ~~ mho/meter 


and T = OS Ea 


exists which is influenced by the pumping theory and experiment. In 


mho/meter. Then ¢, = 0.6. A final possibility 


Chapter V it is shown that the maximum pumping velocities attainable 
for a constant voltage occur at a frequency of 1.5 cps. In the theory 
of Chapter II the frequency dependence of velocity predicts that the 
maximum velocity is strongly influenced by the expression wel (0 + 

(we) ), Which has a maximum when c= Wee For a frequency of 1.5 cps.. 


ano who/meter which occurs at 


= 0.66. At that value of a o =11.2 x 107° mho/meter. For calcu- 


this equation predicts that Na Peston ce LO 


a 


ations in this thesis, these latter values will be used. The other 


possibilities are presented to show the range of uncertainty which mght 


BAG 





be expected as a result of the approximations made in the theory 
of Chapter II. 


te Viscosity 


The viscosity of AROCLOR was measured using three Cannon-Fenske 
viscometers, ASTM Nos. 50, 200 and 400 to cover the temperature range 
of interest. The Cannon-Fenske viscometer measures v, the kinematic 
viscosity. Measurements were made over a range of temperature from 
5° to 68.75 Cs The kinematic viscosity was converted to the absol- 
uve Viscosity, um, by the formula a = pv where 0 is the density. The 
density of AROCLOR was taken from the manufactures ae The results 
of the viscosity measurements are presented in Fig. 5. As with the 
conductivity measurements described earlier, the temperature was related 
to the normalized depth in the channel and the data of Fig. 5 was 
approximated by the analytic function p = p67 where H = Def X 1073 
kg./m.-sec. This seems to be a very good approximation for the 
viscosity. 


De Dielectric Constant 

As mentioned before the dielectric constant for AROCLOR was taken 
from the manufacturers data. The value used was 5.7 measured at 1000 cps. 
"mest for Kinematic Viscosity", (ASTM D445-IP 71) 


X% 
Op. Cit., Monsanto Technical Bulletin No. PL-306, pg. 41 
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Figure 5 Results of measurement of absolute viscosity of 


AROCLOR as a function of temperature and normal- 
ized depth in pumping channel. 
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and 25°C, Since the pumping experiment was run at a mean temperature 
of O°C and 1.5 cps. this value is only a approximation to that which 
existed in the pump. It is probable that the dielectric constant is 
also a function of frequency and temperature, however, this dependence 
was not measured and ¢ was assumed to be constant in the pumping ex- 
periment. 

E onclusion. 

The values of the properties of AROCLOR used for calculation in this 
thesis are shown in Table l. While the viscosity measurements were quite 
straight-forward and gave trustworthy results, the values for o and . are 
more open to question. Further experimentation with more careful measure- 
ment is indicated. There seem to be enough mysteries in the determination 
of the conductivity to justify a considerable expenditure of effort in 


careful experimentation. 


ade 





PROPERTIES OE AROCLOR 


absolute viscosity | 5.7 x 10"kg/m-sec 


= 4.75 x 1077 “inno /meter 


conductivity ise oc 107) imho /meter 


permittivity ge = constant 





Table 1 Properties of AROCLOR used in the evaluating of 
experimental results. 
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CHAPTER IV 


THE EXPERIMENT 


A. Introduction 


The purpose of the experiment performed for this thesis was to 
demonstrate that slightly conducting liquids can be pumped accord- 
ing to the method purposed in Chapter I. Measurements were taken 
to determine how the velocity of the Liquid in the flow depends on 
the voltage of the applied traveling potential wave, and how the 


velocity depends on the frequency of the applied traveling wave. 


Dr Apparatus 


The apparatus used in the experiment is shown in Fig. 6. It 
is a modification of an apparatus used in an earlier experiment in 
EHD pumping of a liquid with a free surface. The liquid to be pumped 
was contained in a re-entrant, circular channel with insulating walls. 
made of plexiglass and a highly conducting bottom made of aluminum. 
The liquid was bound at the top of the channel by an arrangement of 
electrodes to which the traveling potential wave was applied. The 


dimensions of the channel were: 


1) mean circumference 88.6 cm. 
2) depth sen: 
3) width gel cm. 
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Figure 6 Apparatus used in the induction pumping experiment, showing 
traveling wave generator, pumping channel, and hot and cold 
baths. 
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The 40 electrodes which formed the top of the channel were 
individually connected to 40 contacts on a commutator ring. The 
40 contacts were connected in a closed loop with a 22 megohm resistor 
connecting each two adjacent contacts. The commutator was a rotating 
plexiglass bar with a contact mounted at each end. One contact was 
connected through a slip ring to a potential of aa volts and the 
other was connected through another slip ring to a potential of 
EN volts. The commutator was rotated with an angular velocity 
w. In this way, if the commutator was located at some particular 
angular position,ideally the potential on the electrodes varied 
linearly in small steps from Ds at one position at the top surface 
of the channel to =; at the diametrically opposed position. Then 
as the commutator rotated the potential distribution moved around 
the channel as a traveling potential wave with a sawtooth shape. 
The actual wave which was created is shown in Fig. 7. The tendency 
of the wave shape to be more sharply peaked than an ideal sawtooth 
was caused by the loading effect on the power supplies of the liquid 
in the channel. 

The liquid in the channel was AROCLOR 1232. Its properties have 
been described in Chapter III. The experiment requires that a gradient 


in electrical conductivity be established between the top and the bottom 


ne 





ct 


Figure 7 Approximate shape of voltage waveform produced by 
traveling wave generator. Fundemental component 
has a value of 0.66 Hs 
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of the channel. To produce the required gradient advantage was 
taken of the fact that the electrical conductivity of AROCLOR is 

a function of temperature. Consequently, a temperature gradient 
was established in the liquid, between the top and the bottom of 
the channel. The bottom of the channel was put in contact with 

a bath of dry ice in alcohol which kept the bottom of the channel 
at a temperature of approximately =o Or The region above the 
electrodes on the top surface was heated by circulating hot mineral 
oil in a channel above the electrodes. The mineral oil was physically 
Separated from the AROCLOR by the plexiglass ring which also served 
as a mount for the electrodes. Heat from the circulating mineral 
oil was conducted to the AROCLOR by the electrode mounting bolts. 
The mineral oil was maintained at a temperature of approximately 
162°C. By this method, a temperature difference of about 114°C 
was maintained between the top and the bottom of the channel. It 
should be noted, however, that the temperatures measured were the 
temperatures of the alcohol bath and the mineral oil and not the 
temperature of the AROCLOR at the top and bottom of the channel. 
While the heat conduction through the aluminum plate at the bottom 
was probably quite good, it is likely that heat conduction from the 
mineral oil to the AROCLOR was relatively poor. A reasonable 


assumption seems to be that the temperature at the top of the channel 
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was about +50°C and at the bottom, -50°C. This gives a total 
temperature difference of 100% across the channel, or a temp- 
erature gradient of 25°C/cem. , assuming that the variation in 
temperature from the bottom (C = O) to the top (C = 1) of the 
channel is linear. The measurement of the velocity of the fluid 
flow was done by suspending a few small particles of bakelite in 
the AROCLOR and measuring the time that a particle took to pro- 
ceed around the channei. The particles stayed suspended in the 
liquid because a ' - density ‘gradient. was caused by the gradient 
in temperature. The specific gravity of the bakelite particles 

was such that they would sit in the density gradient of the liquid. 
The progress of the particles around the channel was observed 
through the clear plexiglass walls. The particles were about 1/8 in. 
in diameter so that they were massive enough that even if they were 
themselves slightly charged, they would still respond to the fluid 
motion and not more independently of the fluid. 


Ge Procedure 


Two interesting experiments were performed with the apparatus describ- 
ed above. The experiment to measure the dependence of the flow velocity 
on the voltage of the applied potential wave was done at a frequency of 
1.5 cps. for a range of voltages between 7.0 kilovlots and 13.0 kilovolts 
in increments of .5 kilovolts. At each voltage the velocity of a particle 


in the flow was measured several times. At voltage levels below 7.0 
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kilovolts the particle velocities were quite small and somewhat 
difficult to measure. Since the electric field in addition to 
inducing a flow is also tending to leviate the liquid against the 
gravitational force. Above 13.0 kilovolts the liquid became so 
unstable internally that the flow was wiped out by the violent 
local motions of the liquid. The very unstable motion of the 
liquid was easily observable as swirls of liquid of slightly 
different refractive indices which could be seen in the presence of 
a strong light. This experiment was conducted at a frequency of 
1.5 cps, since at a constant voltage, that frequency resulted in 
the largest velocities. Some variations in the velocities measured 
for a given set of conditions were observed and were caused by 
particles sitting at different places in the flow in both the 
WemuLcalana tic horizontal directicons, 

Repeatibility of data was difficult to attain unless the 
AROCLOR was allowed to remain in the channel for several hours 
before the experiment was performed. This may have been due to 
the fact that the AROCLOR was slightly corrosive to the plexiglass 
walls of the channel and some time was necessary to allow this pro- 
cess to achieve some equilibrium. 

Another experiment was performed to find the dependence of 
velocity on the frequency of the applied traveling potential wave. 
This experiment was performed at a constant voltage of 10 kilovolts. 


This voltage level was chosen because the particle velocity was 
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relatively fast and allowed measurements to be made reasonably 
quickly but the liquid was still quite stable with regard to 

violent internal motions. The velocity of the flow was measured 

WA Ag TONG of frequencies from about O cps. 66 5.5 Cps im inver= 
vals of about 0.5 cps. This range of frequency was used because 

of physical limitations imposed by the design of the commutator con- 
tacts. 

In both of the experiments described here, the temperature 
gradient which was established across the channel was constant, 
so that the gradient in electrical conductivity was the same. 
Another experiment which would vary the gradient in conductivity 
could be envisioned but was not performed. The results of the 
experiments measuring velocity as a function of frequency and of 


voltage are shown in Figs. 8 and 9 in Chapter V. 


2372 





CHAPTER V 


RESULTS 


A. Introduction 

In this chapter the results of the experiment described in Chapter 
IV will be presented and discussed in relation to the predictions of 
the theory presented in Chapter II. A program for further experimenta- 
tion will be introduced as well as some ideas for examining the theory 
pertinent to this area of EHD induction in the context of modifications 
to the experiments presented in this thesis. 
B. Experimental Results 

Figure 8 presents the results of the experiment which measured the 
velocity as a function of frequency. The data indicates that for O 
frequency the flow velocity is O. This result is'predicted by the 
theory. The data indicates that maximum vy is achieved atl a ivequenc,, 
SiS cps- and that u gradually decreases at higher frequency. The 
theory indicates that y is strongly dependent on the factor 18/(1+8) 
which may be rewritten as o oK + on Notice that there are also 
meequency terms in the expression Tor aa on frequency in this experiment 


2 
Will probably not, however, be as strong as the dependence on o we/(o_ + 


a 
2 
(we) ). A theoretical curve is plotted on the data in Fig. 8 showing the 


degree of similarity. The exact theoretical prediction of the frequency 
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velocity Vz in cmysec, 


1,4 FLOW VELOCITY 


VS 


TRAVELING WAVE FREQUENCY 
1.2 


1.0 


0.8 


0.6 





connif. X UJG 
Theoretical curve of the form 05 + (wer 





o 1.0 20 3.0 4.0 5.0 frequency 
in C.p.s. 


Figure 8. Results of experimental measurement of flow velocity as 
a function of frequency for N = 10.0 kilovolts with a 
100°C AT between top and bottom of pumping channel. 
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dependence involves a computational problem of modest complexity which 
could be done easily on a digital computer but was not done for this 
thesis. 

As mentioned in Chapter III, the existence of the peak in the 
data of Fig. 8 was used as a key to determining a likely value for 
oo: It is of interest that the value of on predicted by this peak 
is within the range of values generated during the independent measure- 
ment of o by the method of Chapter III. 

Fig. 9 presents the results of the experiment which measured the 
velocity as a function of the voltage. These results are plotted against 
as Ve. The data clearly justifies the theoretical prediction that 
the velocity is linearly dependent on the voltage squared. The con- 
stant of proportionality which is the slope of the estimated line 
through the data points is 1.2 x O E sec/meter. This data 
could be plotted not against the peak voltage of the wave form shown 
in Fig. 7 . of Chapter IV, but instead against the peak value of the 
first fourier component of the voltage. The first fourier component 
has a peak of 0.600 In this case the slope of the line of v, versus 


lO 


(0.66v, Y would be 2.75 x 10° volt”sec /meter. 


Computed for the same wave shape by J.R. Melcher in "Traveling Wave 
Induced Electroconvection", to be published in Physics of Fluids. 
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Figure 9. 
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Results of experiment measurement of flow velocity as a function of 


traveling wave potential V_. 
Measurements made at a frequency of 1.5 c.p.s. 
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Plotted as velocity vs. voltage squared. 





Turning now to the theory, the predicted slope of the line of 
pe versus v may be computed. Using the values for Jo and Sy of 
Chapter III, the value of Y is 2.36. The wave number k for this 
experiment is computed from k = 21r/1 where A, the wave length of 
the traveling wave, is the mean circumference of the channel. By 
this means k is determined to be 7.1 es Since the data was 
taken at a frequency of 1.5 cps., which is the frequency for maximum 
v,» S = waja = 1. With these values and the values for p, and ¢ 
WONO Dhar LIL, Eau. (33) may be numerically integrated. The value 
of ¢ used for the integration is 0.7 which corresponds to the approxi- 
mate position at which the particles used to measure the velocity were 
observed to float. Incidentally, some error is introduced here since 
particles often floated at varying depths but mostly with in a range 
of € = 0.5 to € = 0.8. In any case, performing the numerical integration 
for ¢ = 0.7 by Simpson's Rule gives that the slope of the line for vat 
€ = .7 versus Y” should be -1.17 x D Y . Notice the minus sign. 
After all, the theory predicted that with a positive conductivity gradient, 
inducing @ positive charge in the liquid, the flow should be in the 


direction opposite to the direction of travel of the traveling potential 


wave. The magnitude of the slope certainly agrees closely with the ex- 
beriment. However, in the experiment the liquid proceeded in the wrong 
S 

Ee 





direction, which at first glance is somewhat alarming. 

the resolution of this discrepancy between the theory and the 
experiment may be found in a closer examination of the exact thermal 
conditions in the channel. Unfortunately, in the experiment as per- 
formed and described in Chapter IV, the measurement of the actual tem- 
perature profile in the liquid was never done. The thermal relaxation 
time for AROCLOR is approximately 2 x 10° seconds for distance on the 
order of 3 or 4 mm., which would be a reasonable dimension for the kinds 
of swirling motion observed in the liquid during various experiments. 
This long thermal relaxation time means that when the liquid in a tem- 
perature gradient is stirred, the gradient will not re-establish itself 
in the liquid for some considerable time. As mentioned in Chapter II, 
the electric field imposed on the liquid has a strong tendency to ] evitate 
the liquid against the force of gravity and make it internally unstable. ' 
It may be that the field induces charge in the liquid as expected by the 
positive gradient ino and then the electric forces on this induca charge 
cause it to displace some of the warmer liquid near the top of the channel 
with cooler liquid from the bottom. Then, because of the long thermal 
relaxation time the original temperature gradient is not re-established. 
This is the same as saying that the internal instability of the liquid 
allows the formation of a temperature inversion layer. In this case, a 


region exists in the flow which has a negative conductivity gradient and 
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the induced motion of the liquid is in the positive direction, as the 
theory predicts. The liquid very near the top would still have a positive 
gradient in conductivity as well as that near the bottom. In the experi- 
ment performed, for the most part velocity measurements were only made 
in a range near the mid-depth of the channel. Besides, the very cold 
liquid near the bottom is quite viscous so that velocities there would 
be small in any case. 

It is of great interest that the following observations were made 
during the experimentation. Early in the experiment process, before 
the theory had been fully developed, it was observed that in the pre- 
sence of a gradient of about 30°C, where the liquid was only cooled on 
the bottom but was near room temperature at the top, the liquid was 
observed to flow very slowly in the direction opposite to that of the 
traveling wave. The velocities were small and it was thought that a 
larger temperature gradient would produce greater velocities. A larger 
gradient was applied and indeed larger velocities resulted and in the 
same direction as the traveling wave. Unfortunately, the author's in- 
tuition that the induced flow ought to be in the direction of the travel- 
ing wave, led him to continue working with the larger temperature gradient. 
No meaningful data was taken for the backward flow case. 

During the measurement of pa versus Frequency some particles, which 


were floating just below the mid-depth of the channel, were observed to move 


E 





very slowly backward at the higher frequencies. This is another 
indication of the worth of the theory but this situation has not 


yet been pursued in detail. 


C. Further Work 

This experiment has shown that electroconvection can be induced. 
A theory now exists which correlates with the experiment in some ways. 
Barcher experiments are indicated, however: Ii the temperature eradı nv 
which actually exists in the Liquid under the conditions described in 
Chapter IV were carefully measured, the results of this experiment 
could be better interpreted in relation to the theory. 

Another experiment might be done to demonstrate the flow in the 
backward direction. The distribution of positive charge which results 
im the Liquid with a positive conductivity gradient is different than 
the distribution of negative charge which results with a negative con- 
ductivity gradient. This may be seen in Equ. (34) of Chapter II. 
Therefore, if the sign of the conductivity gradient changes but the 
absolute magnitude remains the same, both the sign and the magnitude of 
the velocity will change. That is pumping backward is not exactly the 
Kovero of pumping Torward. 

The most interesting extension of the theory which should be made 
is to do what might be called a multi-region analysis. This would only 


be a matter of considering the liquid in the channel sliced into different 
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regions. The conductivity approximation o = om + VEG) could 
De very econ beer e cion WIDA proper sized slices oF would 
be sufficiently larger than g, so the approximation in Equ.(9 ) that 


09, but V o = o, would be very good. The appropriate boundary con- 


1 
ditions would be to match fields and velocities at each interface be- 
tween regions. This problem would require computer solutions. 

The theory of Chapter LI is the theory for the long wave limit 
which is the case for which the term 2dk << 1/1 + joa Ine short wave 
limit which is the reverse inequality is presented in a very brief form 
in Appendix B. An experiment might also be performed where this limit 
Aad to test she Theser tCAToOuen ICS entire ocope. 

Another major area of interest was mentioned briefly in Chapter III. 
The schemes discussed in this thesis for EHD induction pumping have possi- 
bilities as methods for further definition of the mechanism for conduction 


in liguids and for the investigation of electrically induced internal 


instabilities in liquids. 


nice 





REFERENCES 


Chandrasekhar, S., Hydrodynamics and Hydromagnetic Stability, 
London 1901. 


Pano R M Cnu, Leer and Adier MR: B., Electromagnetic 
Fields Energy and Forces, New York, 1960. 


Melcher, J. R., "Traveling Wave Induced Electroconvection”, 
SO be published 1m Physics of Pude. 


Woodson and Melcher, Class notes for Fields, Forces, and 
Motions., M.I.T., Cambridge, Massachusetts, 


1964. 


Monsanto Chemical Company, "Aroclor Plasticizers", Technical 
Bulletin No. PL - 306, St. Louis, 1960. 


1965 Book of ASTM Standards, Part 17, Section (ASTM D 445 - 
Terran, Test for Kinematic Viscosity- 


nTa 








APPENDIX A 


SOLUTION FOR A SINGULAR CONDUCTIVITY GRADIENT 


A problem similar to the bulk EHD induction pump 
would be a pump operating with two liquids of differing 
o's and e's. Consider a device with the same electrical 
arrangement as that described in Chapter IV but with a 
liquid in the upper region of parameters O, and E, and 
a liquid in the lower region of o, and €}. (See Figure 
10). In this case the induced charges will relax to the 
interface forming a sheet of charge at y = 0. Solving 
Maxwell S equations for the electric case will result 


in potential s in the two regions of the form: 


ð = Re [A sinh ky + B cosh ky | eot- kx) (1) 
the resulting field will be 


A 
E 
y 


-kj A cosh ky + B sinh ky | (2) 


A 
E 
X 


jk [A Sina ky + B cosh ky] (3) 
The boundary conditions are: 


aty sa b =N 





j (wt -kx) 
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Figure 10. Coordinate system for the two fluid problem with 
@ Singularity in conductivity gradient at the 
interface. 
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at y = 0 Ex, = EX and 


S 3 
EY T 9EY; t 5E fe Fy, - &pFy,) t U SF Ce Ey -e Fy) = 0 
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where the u and £ indicate the upper and lower regions 
and U is the velocity in the x direction at the interface. 
Matching these conditions with the expressions for the 
potential and the fields in each region, equations (1), 


(2), and (3) gives: 


We A sinh ka + B cosh ka (4) 
O=B -B (5) 


0 = (o, + jw-kU)E,) AL - Coti lo-kU)e,) A, (6) 


O=- A sinh kb + B, cosh kb (7) 


The time-average electric traction at the surface will 


be given by the Electric Stress Tensor and will be: 


T eee ei 8 Ey e 8 
< Xy gS ea on ee ee > 


=90= 





After solving for the A's and B's from equations (4) to 
(7), then combining equations (2) and (3) with (8), the 


time-average traction becomes: 


Sn kb sinh kb)e (w-kU)(o,€ -o €,) 
< T E O u Lu ou 8 


ay 2 [0 Cosh ka sini kb f o, sinh ka cosh Kb |? 


£ 


T C | € sinh ka cosh kb + E, cosh ka sinh kb | ¿ 


(9) 


If equation (9) is taken in the limit where oe O and 


E then the result is the same as that derived by 


Ye 
Melcher for the single liquid with an air interface. 


That is, for w > kU the traction is positive. If, however, 


€ = Ey but oi > y conditions which would approximate 


4 5 
those under which F experiment of this thesis was done, 
then the time-average electric traction is negative. In 
that case, U would surely be negative. The sign of the 
charge induced on the interface would be negative. Thus 
the result derived here for an exact solution agrees with 
the solution for the continuous gradient in conductivity 


of Chapter II. 


*Melcher, op. cit. 





APPENDIX B 


Solution For Small Conductivity Gradient 


From Equ. (28) of Chapter II a general expression for the time- 
average electric traction may be written in terms of the a's and b's, 


before taking thelimit that 2dk<< 7/1+jS. That expression is: 


2 E e r = 
ekV e = 1) b sin b.C-b.Sinh b € + a. (cosh pC idos D.C) 
<T E O T Ta T al e u 
xy La (cosh by = cos b.) 
r i 
If instead of the above limit one takes the limit that 2dk>> 7/1+js 


then the a's and b's are (from Equ. (26) ): 
2 
a, =b, = 03  b, =2dx, O (1 + S°) 


Substituting these values into the above expression for ae gives: 


E 2 
> = — | = 
Ty ekV_ 18 (cosm 2dke =) 


Ya es (cosh 2dk -1) 


Now if 2dk<<. this may be further approximated as: 


2 
E vé TSG 
E mn 
ha (1+87) 


| E an 
With this new TS and p = pS the velocity Equ.(30) may be integrated 


Dr 





to give: 


Tí 2 
a e k v 18 6 de c.a a e 
xX we AAA 


Wa. eee 7 u 8 
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When C- and C, are evaluated for the no-slip condition at G = O and 


( = 1 the result is; 


Dee) 


aan 2 
= 281, (1 +S) 





This equation for the limit of small conductivity gradient shows that 
in the limit where S becomes very large that the velocity — O which 
means that the relaxation time is Goo long for The device To work at 
high frequencies. The equation also predicts that vy, <0 forn- O, 


and vy =O Lore] < O- 


-53= 
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